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Abstract 

The solutions of the Wigner-transformed time-dependent Hartree-Fock-Bogohubov 
equations are studied in the constant-A approximation. This approximation is 
known to violate particle-number conservation. As a consequence, the density fluc- 
tuation and the longitudinal response function given by this approximation contain 
spurious contributions. A simple prescription for restoring both local and global 
particle-number conservation is proposed. Explicit expressions for the eigenfrequen- 
cies of the correlated systems and for the density response function are derived and 
it is shown that the semiclassical analogous of the quantum single-particle spectrum 
has an excitation gap of 2A, in agreement with the quantum result. The collective 
response is studied for a simplified form of the residual interaction. 
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1 Introduction 



The problem of extending the Vlasov equation to systems in which pairing 
correlations play an important role has been tackled some time ago by Di 
Toro and Kolomietz [1] in a nuclear physics context and, more recently, by 
Urban and Schuck [2] for trapped fermion droplets. These last authors derived 
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the TDHFB equations for the Wigner transform of the normal density matrix 
p and of the pair correlation function k (plus their time-reversal conjugates) 
and used them to study the dynamics of a spin-saturated trapped Fermi gas. 
In the time-dependent theory one obtains a system of four coupled differential 
equations for p, and their conjugates [2] and, if one wants an analytical so- 
lution, some approximation must be introduced. Here we try to find a solution 
of the equations of motion derived by Urban and Schuck in the approximation 
in which the pairing field A(r, p, t) is treated as a constant. It is well known 
that such an approximation violates both particle-number-conservation and 
gauge invariance (see e.g. sect. 8-5 of [3] and [4]), nonetheless we study it 
because of its simplicity, with the aim of correcting the final results for its 
shortcomings. Moreover, the constant-A approximation is not satisfactory for 
describing long wavelength pairing modes in a large system. Such modes have 
frequencies which are much less than the pairing frequency A//i and for their 
study it is essential to use a self consistent theory where the gap A is related 
to the pair density k through the pairing interaction. The phases of A and 
K are particularly important because they describe the coHective superfiuid 
currents. On the other hand nuclei are small systems. Shell gaps are large 
compared with A, or equivalent ly giant resonance frequencies are large com- 
pared with the pairing frequency. The constant-A approximation is much more 
reasonable in such systems. 



In Sect. 2, the basic equations are recalled and reformulated in terms of the 
even and odd components of the normal density P- In Sect. 3, the static limit 
is studied by following the approach of [5] and the constant-A approximation 
is introduced. In Sect. 4, the simplified dynamic equations resulting from the 
constant-A approxmation are derived and their solutions are determined in lin- 
ear approximation. In Sect. 5, these solutions are studied in a one-dimensional 
model and the problem of particle-number conservation is examined in detail. 
By studying the energy- weighted sum rule (in the Appendix), we find that the 
constant-A approximation introduces some spurious strength into the density 
response of the system. A simple prescription, based on the continuity equa- 
tion, is proposed in order to eliminate the spurious strength. The resulting 
strength function gives the same energy- weighted sum rule as for the uncor- 
related systems. In Sect. 6, the general solution found in Sect. 4 is re- written 
for spherical systems, where the angular integrations can be performed ex- 
plicitly, leading to expressions containing only radial integrations. In Sect. 
7, the collective response function of spherical nuclei is derived for a simple 
multipole-multipole residual interaction. In Sect. 8, the quadrupole and oc- 
tupole channels, that are the ones most affected by the pairing correlations 
are shown explicitly. Finally, in Sect. 9 conclusions are drawn. 
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2 Basic equations 



We assume that our system is saturated both in spin and isospin space and 
do not distinguish between neutrons and protons, so we can use directly the 
equations of motion of Urban and Schuck. 

We start from the equations of motion derived in Ref. [2] for the Wigner- 
transformed density matrices p — p(r, p, t) and k — K{r, p, t), with the warning 
that the sign of k, that we are using agrees with that of Ref. [1], hence it is 
opposite to that of [2]. Moreover we find convenient to use the odd and even 
combinations of the normal density introduced in [2]: 



Pev = ^[p{r,p,t)+ p{r,-p,t)], (1) 
Pod = ^ [p(r , P, t) - p{r, -p, t)]. (2) 



Thus, the equations of motion given by Eqs.(15a...d) of Ref. [2] read 

ihdtPev^ifi'{h,Pod} - 2ilm[A*{r,p,t)K] (3) 
ihdtPod = 'ifi{h,Pev} + ihRe{A*{r,p,t),K,} (4) 
ihdtK = 2{h - p)k - A(r, p, t){2pev - 1) + i^{A(r, p, t),Pod} ■ (5) 

Here h is the Wigner-transformed Hartree-Fock hamiltonian /i(r, p, t), while 
A(r, p. t) is the Wigner-transformed pairing field. Since the time-dependent 
part of K, is complex, k, = K^+iKi, the last equation gives two separate equations 
for the real and imaginary parts of k. 

Moreover, from the supplementary normalization condition ([6], p. 252) 

= n (6) 

satisfied by the generalized density matrix TZ, the two following independent 
equations are obtained: 



Podl^ + i^{Pev,K}^0, (7) 
PeviPev - 1) + Pld + KK*^0. (8) 

We shall use the equations of motion (3-5), together with these equations, as 
our starting point, but first we notice that, in the limit of no pairing, both A 
and K vanish, the third equation of motion reduces to a trivial identity, while 
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the first two give the Vlasov equation for normal systems, expressed in terms 
of the even and odd components of p: 



dtPev^{h,Pod} , (9) 
dtPod^{h,Pev} ■ (10) 

A solution of the linearized Vlasov equation for normal systems (i. e. without 
pairing) has been obtained in Ref. [7] and our aim here is to study the changes 
introduced by the pairing interaction in the solution of [7]. 

Moreover, before studying the time-dependent problem, it is useful to look at 
the static limit. 



3 Static limit 



In this section we follow the approach of Ref. [5] . At equilibrium we have 



Pev = Po{r,p), (11) 

Pod = 0, (12) 

/i = /io(r,p), (13) 

K = Ko{r,p) (14) 

Ao = Ao(r,p) (15) 

and equations (3-5) give 

= -2ilm(A*Ko), (16) 

= ih{ho,po} + ihRe{A*,Ko}, (17) 

= 2(/io-/i)«o- Ao(2po-l), (18) 
while Eqs. (7, 8) give 

t^{po,f^o} = 0, (19) 

Po(po-l) + |«:o|' = 0. (20) 



Equation (16) is satisfied if we assume that Aq and are real quantities, 
while Eqs. (18) and (20), taken as a system, have the solution: [5] 

Po(r,p)^-(l- ^^^^p^ ) (21) 
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I ^ Ao(r,p) 

with the quasiparticle energy 

E{v, p) = V^2(r,p) + (/,,(r,p)-;,)2 . (23) 

It can be easily checked that Eqs. (21, 22) satisfy also Eqs. (17) and (19), that 
is 

{/io,Po} + {Ao,«:o} = (24) 

and 

{Po,«o} = 0. (25) 



The (semi) classical equilibrium phase-space distribution is closely related to 
Po(r,p): 

4 

/o(r,P) = ^^^Po(r,p) (26) 



and the statistical factor 4 takes into account the fact that there are two kinds 
of fermions. 

The parametrer is determined by the condition 

A = ^ rfrrfp/o(r,p), (27) 

where A is the number of particles. This integral should keep the same value 
also out of equilibrium (global particle-number conservation) . 

3.1 Constant-^ approximation 



In a fully self-consistent approach, the pairing field A(r, p, t) is related to 
«;(r, p, i), however here we introduce an approximation and replace the pairing 
field of the HFB theory with the phenomenological pairing gap of nuclei, hence 
in all our equations we put 

A(r,p,i) Ao(r,p) a; A = const , (28) 
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with A IMeV. 

In the constant-A approximation the equihbrium distributions become 

%W = -^ (30) 

and the quasiparticle energy 



E{e) = VA2 + {e- , (31) 

with 

e = ho{r,p) = ^ + Vo{r) 
the particle energy in the equihbrium mean field. 

In the following we shall use the relation: 

E^e) dpoje) 

«o(6) = ^ . (32) 

4 Dynamic equations 



Always in the approximation where A is constant and real, the time-dependent 
equations (3-5) become 

ihdtPev = iH^, Pod} - 2iAlm{K) (33) 

ifidtPod = ih{h,Pev} (34) 
ihdtn^2{h- p)n- l^{2p^^-l) . (35) 



This is the simplified set of equations that we want to study here. The sum of 
the first two equations gives an equation that is similar to the Vlasov equation 
of normal systems, only with the extra term — 2iAIm(«;). This extra term 
couples the equation of motion of p with that of k, thus, instead of a single 
differential equation (Vlasov equation) , now we have a system of two coupled 
differential equations (for p and Ki). 

Our aim here is that of determining the effects of pairing on the linear response 
of nuclei, thus we assume that our system is initially at equilibrium, with 
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densities given by Eqs. (29,30), and that at time t — a. weak external field 
of the kind 



6V'^\r,t) = /3S{t)Q{r) (36) 



is applied to it. This simple time-dependence is sufficient to determine the 
linear response of the system. In a self-consistent approach, we should take 
into account also the changes of the mean field surrounding each particle 
induced by the external force and consider a perturbing hamiltonian of the 
kind 

Sh = SV^' + , (37) 



however we start with the zero-order approximation 

5h = 5V^^^* (38) 

and will consider collective effects in a second stage. 

Since we want to solve Eqs. (33-35) in linear approximation, we consider small 
fluctuations of the time-dependent quantities about their equilibrium values 

and neglect terms that are of second order in the fluctuations. Hence, in Eqs. 
(33-35) we put: 

h^ho + Sh, (39) 

Pev^Po + SPev: (40) 

Pod^SPod, (41) 

K — Kq + 5k — Kq + SKr + iSKi . (42) 

Then, the linearized form of Eqs. (33-35) is 

ihdtSpev = ih{ho, Spod} - 2iA5Ki (43) 

ihdtSpod = ifi{ho, Sp^v} + ih{Sh, po} (44) 

—hdtSKi — 2(e — p)SKr + 2KoSh — 2ASpev , (45) 

hdtSur — 2(e — p)5K,i . (46) 

Taking the sum of the first two equations gives 

indtSp{r,p,t)= (47) 
ih{ho, Sp{r, p, t)} + ih{Sh{r, p, t),po} - 2iASni{r, p, t) , 
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which can be regarded as an extension of the hnearized Vlasov equation stud- 
ied in [7]. 

In order to make the comparison with [7] easier, from now on we change the 
normahzation of the phase-space densities and define 



/(r,P,^)=(^P(r,P,t) (48) 
4 

x{r,P,t)^-^^^K{r,p,t), (49) 
moreover, we put 



In terms of the new functions Eqs. (47) and (45) read 



ihdtSf{r, p, t) = ih{ho, Sf{r, p,t)} + zh{6h{r, p, t), fo} 

-2iA5xiir,p,t), (51) 
-hdtSxi{r, p, t) = 2(e - ii)Sxr{r, p, t) + 2xoSh{r, p, t) 

-2A5fe,{r,p,t). (52) 

The function f^^ is given by the obvious extension of Eq. (1). In order to get 
a closed system of equations, we still need an extra equation for 5xj.(r, P, ^)- 
This can be obtained from the linearized form of the supplementary condition 
(8) that reads 

Spev{2po - 1) = -2KodKr , (53) 



or 

5Kr{r, p, t) = ^ ^^^^^h pev{r, p, t) = -^—^5pevir, p, t) . (54) 

ZKq Z\ 

The last expression has been obtained with the help of Eq. (18). In terms of 
the new functions / and the last equation reads 

SXr(r, p, t) = -^-^Sfev(r, p, t) . (55) 



Equations (51, 52) and (55) are the set of coupled equations for the phase- 
space densities that we have to solve. 
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Replacing Eq. (55) into Eq. (52), and using Eq. (30), gives the following system 
of coupled differential equations: 



dtSfir, p, t) = {ho, 5f} + {5h, fo} - 2^Sxi{r, p, t) , (56) 
dtSxiir, p, t) = P, t) + 5f{r, -p, t)] - 2^Sh{r, p, t) . (57) 

Taking the Fourier transform in time, gives 



-iu;5f{r, p, u) = {ho, 5f} + {5h, fo} - 2^5xi{r, p, u) , (58) 

-iu}6xi{r, p, a;)) = -^[^/(r, P, ^) + Sf{r, -p, a;)] 

-2^5h(r,p,u;), (59) 

or (for uj ^ 0) 

-iu5f{v, p, u) + {5/, /lo} = -i^d'^]^[5f{r, p, a;) + 5/(r, -p, a;)] 

+ F'(e) [{5h, ho} + iujdFSh] , (60) 

with 

- f^)' (61) 



and 



Q(e) = 2^ . (62) 

h 



This frequency plays a crucial role in our approach, its minimum value is 

2A/;i. 

In Eq. (60) we have used the relation {jo, ^K} ~ F'{e){ho, Sh} as well as Eq. 
(32). 

By comparing Eq. (60) with the analogous equation for normal systems 

-iu;Sf{r,p,u;) + {Sf, ho} = F'{ho){Sh, ho} , (63) 

we can see that the only effect of pairing in the constant-A approximation is 
that of adding the terms proportional to (P. 
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The normal Vlasov equation (63) can be solved in a very compact way by using 
the method of action-angle variables [8] , [7] . In that approach one expands 

5/i = ^5/i„(I)e--*, (64) 



where I and $ are the action and angle variables, respectively. Moreover 

Sf{r,p,cu) = J26Ul,u)e^--^ (65) 

n 

and 

{5/, ho} = E^(n • c^)5/„(I,cu)e-* , (66) 

n 

where the vector Co has components 

^a=-Kl-- (67) 

Then Eq.(63) gives 

The (zero-order) eigenfrequencies of the (normal) physical system are 

cUn = n ■ a; . (69) 



Here want to use the same method to solve the more complicated equation 
(60). Since that equation contains also the function 5f{r, — p,a;), we need also 
the analogous equation for this other quantity: 



-iu}5f{r, -p, u) - {Sf, /io}r,p = -iu}d^^[Sf{r, p, tj) + Sf{r, -p, u)] 

+ F'{e) [-{5h, /io}r,p + ioJ(f5h] . (70) 

By expanding (5/(r, p, a;) and 5/(r, — p, a;) as 

5/(r,±p,u;)=5:<5/±(I,u;)e-*, (71) 
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Eqs. (60) and (70) give 



[-^(1 --)+ + LU-Sf- = F'(6)K + ujd']dh^ , (72) 

[^y]5/n + [-^(1 - y) - ^n]^/„ = F'(e)[-a;„ + a;ci2]5/i„ , (73) 

which is a system of two coupled algebraic equations for the coefficients 5/+ 
and Sf~. Its solution is 



Sf: = ^|^F'(6)5/.„ , (74) 

5/„ = ^^^f^F'{e)Sh^ , (75) 
where 

^^ = ^n + ^^'(e) (76) 



are the (squared) eigenfrequencies of the correlated system. These eigenfre- 
quencies are in agreement with the enegy spectrum of a superfiuid infinite ho- 
mogeneous Fermi gas (see e. g. Sect. 39 of [10]) and they lead to a low-energy 
gap of 2A in the excitation spectrum of the correlated systems. However, 
as anticipated, we expect problems with particle-number conservation. These 
problems are better discussed in one dimension, where formulae are simpler. 



5 One-dimensional systems and peirticle-number conservation 



In one dimension, Eq.(60) reads 



-iuj5f{x, p, uj) + xdjf - 



dx 



dJf 



(77) 



-iujd^h6f{x,p,uj) + 6f{x, -p,uj)] + F'{e){^dJh + iujd^Sh) . 
2 m 



In zero-order approximation 5h{x,p,u!) — (3Q{x), moreover, in one dimension 
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The vectors u and $ have only one component: 

^0(6) = ^ (79) 



and 

^x) = uJor{x) , (80) 

with 



1 

r{x) = / rfx'— — , (81) 



Xl 



v{e,x) = ^^[e-Vo{x)]. (82) 

The time T(e) is the period of the bound motion of particles with enegy e 
in the equilibrium potential well Vo{x): T = 2t{x2). The points Xi^2 are the 
classical turning points for the same particles. Instead of the action variable 
/(e) = ^ § dxp{e,x), it is more convenient to use the particle energy e as 
constant of motion. As pointed out in [7], the range of values of r can be 
extended to the whole interval (0, T), by defining 

T " 



when T > ^. With this extension, the angle variable ^{x) takes values between 
and 27r, as it should. 

In one dimension, Eqs. (74, 75) give 

5f^{e,u;)^^p^F\e)5h^, (84) 



with 

6hn = PQr 



v{e, x) 



^-iuJnT{x) 



fi^jdx 

P^hx-^cosMx)], (85) 
1 J v(e,x) 

Xl 
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The frequencies a;„ are the eigenfrequencies of the uncorrelated system: 

ujn = rujjQ , (86) 

while CUn are the new eigenfrequencies modified by the pairing correlations: 



cUn = ±Jojl + n^{e). (87) 



Note that, since Sh-n = Shn, then 5/„ = 

By using the solutions (84), we can also obtain an expansion for the even and 
odd parts of Sf: 



dfev{x, e,(jj) = ^ An{uj) COSmL!oT{x) , (88) 

n=0 

oo 

Sfod{x, e,uj) = ^ Bn{oj) smnujQT{x) , (89) 



n=l 

with 



A^iu^) = -^F'ie)5K , (90) 



B^{uj) = tuj^^^F'{e)6h'^ (91) 
and 

Sh'„ = 2Shn , n^O, (92) 

= Shn, n=0. (93) 

Note that, while Bn=o{i^) — 0, we have An=o{uj) 7^ 0, and this fact leads to 
an unphysical fluctuation of the number of particles, induced by the applied 
external field. These fiuctuations are given by 

5A{uj) = J dx5Q{x,uj) , (94) 

where 5q{x, uj) is the density fluctuation at point x: 

5q{x,u;)^ J dp5f{x,p,u;) ^2 J 8fev{x,e,u) . (95) 
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Equation (88) gives 



T/2 



SA{lj) = An{Lj) / drcosnouoT. (96) 

n=0 i 



Since the integrals /g drcosnwoT vanish when n 7^ 0, the term with n = 
is the only one contributing to this sum, thus givig an unphysical fluctuation 
of the number of particles. This problem could be solved simply by excluding 
the mode n = from the sum in Eq. (88), however this would not be sufficient 

to solve all problems with particle-number conservation, since we can easily 
check that the solutions (88, 89) do not satisfy the continuity equation 

iuQ{x, u) = 9j,j(x, u) . (97) 



The density fluctuation involves only the even part of 5/, while the current 
density j(x, a;) involves only the odd part: 

3{x,uj) = / dp—Sf(x,p,uj) = 2 / de5fod{x,e,u) . (98) 
J m J 



The fact that the continuity equation is violated is a very serious shortcom- 
ing of the constant-A approximation. However, since we have seen that this 
approximation leads to very simple equations and to rather satisfactory ex- 
pressions for the eigenfrequencies of the correlated systems, we still use it, 
but with the following prescription: when calculating the longitudinal response 
function, the density fluctuations should be evaluated by using Eq. (97), instead 
of Eq, (95). Then, the density fluctuations (95) should be replaced by 

2 r 

Sg{x,Lj) = — dedJfod{x,e,u;) . (99) 

lU J 



In practice we are proposing to evaluate the longitudinal response function 
in terms of the transverse response function. It is well known that also the 
more familiar BCS approximation gives a more accurate description of the 
transverse response (see e.g. sect. 8-5 of [3]). In the Appendix we show that the 
longitudinal response function resulting from the present prescription satisfles 
the same energy- weighted sum rule as the uncorrected response function. This 
would not necessarily happen if, instead of changing only the even part of 5f, 
we had modifled also its odd part. 

It is interesting to see how the solutions (84) are changed by our prescription. 
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By using Eq. (89) ioifod, Eq. (99) gives 

/dc - 
— rSfev{x,e,Uj) , (100) 
V[€, X) 

with 

5fev{x, e,uj) = Yj ^n(<^) cosnuJoT{x) , (101) 



n=0 



and 



A^ico) = ^Br^icu) , (102) 



(note that An=o{uj) = 0). Then 

Sf{x, ±p, u) = Sfev{x, e, a;) ± Sfod{x, e, a;) (103) 



and Eq. (84) is replaced by 

Sme,u;)^4^^F'ie)SK. (104) 



By comparing this expression to Eq. (84), we can see that the fluctuations 
of the phase-space density given by the constant-A approximation contain an 
extra contribution that we identify as spurious: 

Sf^{e,u;) = Sf^{e,u;) + Sfr'{e,u;) , (105) 

with 

Sfr^ie, CO) = ^P^F'{e)5K ■ (106) 

The spurious character of 5/**'"'' is suggested also by sum-rule arguments 
(see Appendix). The term /^'""(e, a;) contributes to all modes of the density 
strength function: the contribution to the mode n — Q gives a fluctuation of 
the particle- number integral (global paticle- number violation) , while the other 
modes give a spurious contribution to the density strength function, increas- 
ing the sum rule and violating the continuity equation (local particle-number 
violation). Note that the spurious contribution (106) affects only the even part 
of the pase-space density, not the odd part. 
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6 Spherical Systems 



The method of action-angle variables gives a very compact solution of the lin- 
earized Vlasov equation both in the uncorrelated and correlated cases, however 
it may be useful to make a connection between the results given by this method 
and the more explicit treatment of spherical nuclei given in [7]. For uncorre- 
lated system this has been done in [11]. Here we follow that approach in order 
to derive useful expressions for correlated spherical systems. The components 
of the vecor n are (ni, 77,2, n^), the first point to notice is that, because of the 
degeneracy associated with any central-force field, the vector u) has only two 
non vanishing components: 



With A we denote the magnitude of the particle angular momentum. According 
to Eq. (69), the eigenfrequencies of the uncorrelated system are [7] 



(107) 



'n3,n2(e, A) = ns^o + n2UJ^ , 



(108) 



while Eq. (76) gives the correlated eigenfrequencies 




(109) 



In three dimensions, the Fourier coefficients analogous to (85) are 




(110) 



The external field Q{y) can be expanded in partial waves as 



Q(r) = E<5L(r)lLM(f), 



(111) 



LM 



givmg 



(112) 



LM 



with [11] 



QiLM)^ E yLN{^,^)d'itNWM,nMn.Qn, 

AT— T ^ 



N=-L 



(113) 
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By using this last equation (and changing ris — > n), the expansion (71) be- 
comes 



oo L L oo 

5/(r,±p,a;) = E E E E (114) 

L=0 M=-L N=-L n=-oo 

Sf^^ie, A, a;)e^^-W (©^(a, P', 7)) V^^ (|, |) , 

with 

<^/."^(e,A,a;) = %^^/5F'(e)g^^ (115) 



and Q^jv the semiclassical hmit of the radial matrix elements: 
1 /■ dr 



2 r dr 



/cLt 
— - cos[0„^(r)]g4r) . (116) 

Here T is the period of radial motion, Vr{r) the radial velocity 



"'W = V|(-W-^) (117) 



and the phases (f)nN{r) are given by 

0niv(r) = a;„ivT(r) - Nj{r) , (118) 



where 



and 



r dr' X 

7(r)=/4^A- (120) 
J Vr{r) mr^ 



The frequencies ujq and o;,^ are given by 
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^--^)- (122) 
The Wigner rotation matrix elements in Eq. (114) are given by [9] 

P^^(a, 7) - e-^"dU(/5')e-^^" , (123) 

where (a,/5',7) are the Euler angles introduced in [7]. 

On the basis of the discussion in Sect. (5), we expect that the solution (114) 
will contain some spurious strength introduced by the constant-A approxima- 
tion. In order to eliminate the spurious contributions, we should replace the 
coefficients (115) with 

S^e^ A, a;) = %^^/3F'(e)g^^ . (124) 

These modified coefficients allow us to obtain the modified zero-order propa- 
gator 

Dl{r, r', uj)= j deF'ie) / A £ ^2 , (125) 



with 



d^Nir,r')= (126) 

(47r)^ .TT TT p / ~2uJnN \ / ^nN \ COS 0nAr(r) COS (j)nN{r') 

2L + l'^'^4' 2'" V T ^^OnN^ rV(r) r'V(r') 
and the corresponding response and strength functions: 

TZliu) = J drdr'r^QLir)Dl{r, r', cuy^QLir') , (127) 
SKlo) = --ImlZKuj) . (128) 

TT 

For multipole response: Qiir) — r^. 

For normal systems, the zero-order propagator D^{r,r',uj) is given by Eqs. 
(125) and(126) where CUnN is replaced by uJnN and F'{e) is proportional to a 
5-function[7]. 
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7 Collective response 



Up to now, we have been concerned only with the zero-order approximation, 
which corresponds to the single-particle approximation of the quantum ap- 
proach. In this approximation, the perturbing part of the hamiltonian is given 
only by the external field, while a more consistent approach would require 
taking into account also the mean-field fluctuation induced by the external 
force, so that 

5h = (5y"^*(r, u) + (5y^"*(r, u) . (129) 



In the Hartree approximation, 

5y^"*(r, u)^ j dr'v{r - r')Sg{r', u) . (130) 

where v{r — r')is the (long-range) interaction between constituents. 
For consistency, we take 

SV'''\r,u;) = J dr'v{r-r')dg{r',uj), (131) 



then the collective propagator for correlated systems satisfies the same kind 
of integral equation as for normal systems [7]: 



DL{r,r',uj)= (132) 
Dl{r,r',u) + j dxx"^ j dyy'^DKr, x,u)vLix,y)DL{y,r' ,cj) . 

Here VL{x,y) is the partial- wave component of the interaction between par- 
ticles. We assume that this interaction can be approximated by a separable 
form of the kind 

VL{x,y) = KLX^y'^ , (133) 



where /t^ is a parameter that determines the strength of the interaction. Then, 
the integral equation (132) gives an algebraic equation for the collective cor- 
related response function 

1Zl{ou) = J drr'^r^ J dr'r'\'^DL{r,r' ,u) 
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leading to the expression 



HlM = , , . (134) 

1 - klTZi{lj) 



8 Results 



Here we compare tlie multipole strength functions given by our simphfied 
model of pairing correlations with that of the corresponding uncorrelated sys- 
tem. This comparison is made for the quadrupole and octupole strength func- 
tions, since these channels are the ones that are most affected. 

The static nuclear mean field is approximated with a spherical cavity of radius 
R — 1.2 A3 fm and the A nucleons are treated on the same footing, i. e., we do 
not distinguish between neutrons and protons. Moreover, we chose A — 208 
for ease of comparison with previous calculations of uncorrelated response 
functions [12,13]. Shell effects are not included in our semiclassical picture and 
the results shown below should be considered as an indication of the qualitative 
effects to be expected in heavy nuclei. For the uncorrelated calculations, the 
Fermi energy is determined by the parametrization chosen for the radius as 
ep ~ 33.33 MeV, while for the correlated case, the parameter fi is determined 
by the condition (27); with the value of A = 1 MeV used here, the value of 
is practically coincident with that of ep, so we have used /i — ep — 33.33 MeV 
in the calculations below. Moreover, the small parameter e appearing in Eq. 
(125) has been given the value £ = 0.1 MeV. This value is chosen to simplify 
the evaluation of the response function by smoothing out discontinuities in its 
dependence on lu. 

In the evaluation of the collective response, the value of parameters kl is 
the same as in [12,13], that is: K2 — —1 x 10~^ MeV/fm^ and K3 — —2 x 
10-5 MeV/fm^. 



8.1 Quadrupole response 



Figure 1 show the longitudinal quadrupole strength function evaluated in the 
zero-order approximation (corresponding to the quantum single-particle ap- 
proximation). The dashed curve shows the uncorrelated response evaluated 
according to the theory of [7], while the full curve shows the result of the 
present correlated calculation. As we can see the effect of pairing correlations 
on this zero-order strength function is rather small, however, since pairing af- 
fects also the real part of the zero-order response function, in Fig. 2 we plot 
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Fig. 1. Quadrupolc stcngth function in zero-order approximation. The dashed curve 
gives the response of a normal system ol A = 208 nucleons contained in a spherical 
cavity, while the solid curve includes the effects of pairing correlations in constant-A 
approximation. 
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Fig. 2. Collective quadrupole strength function showing the giant quadrupole res- 
onance. The solid curve involves also pair correlations, the dashed curve has no 
pairing. 

also the collective strength function given by Eq. (134). Again, the effect of 
pairing correlations is very small, in agreement with the results of [1]. 
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Fig. 3. Same as Fig.l, at low excitation energy. 

The main difference between tlie uncorrelated and correlated responses occurs 
at small excitation energy, Fig. 3 sfiows a detail of Fig. 1 at low excitation 
energy. The correlated strength function displays a gap of about 2 A , the 
very small strength extending below 2 MeV is entirely due to the finite value 
of the small parameter e used in the numerical evaluation of the propagator 
(125). 



8.2 Octupole response 



Figures 4 and 5 show the zero-order and collective octupole strength functions, 
both correlated and uncorrelated. As we can see, in this case too the effect is 
rather small. 
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Fig. 4. Zero-order octupolc strength function. The sohd curve involves pair correla- 
tions, the dashed curve has no pairing. 
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Fig. 5. Collective octupole strength function. The solid curve involves pair correla- 
tions, the dashed curve has no pairing. 

9 Conclusions 

The solutions of the semiclassical time-dependent Hartree-Fock-Bogoliubov 
equations have been studied in a simplified model in which the pairing field 
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A(r, p, t) is treated as a constant phenomenological parameter. Such an ap- 
proximation is known to violate some important constraints, like global (particle- 
number integral) and local (continuity equation) particle-number conserva- 
tion. In a linearized approach , we have shown that the global particle-number 
violation is related only to one particular mode of the density fluctuations, 
while the violation of the continuity equation gives a spurious contribution 
to all modes of the density response. Both global and local particle-number 
conservation can be restored by introducing a new density fluctuation that is 
related to the current density by the continuity equation. This prescription 
changes the strength associated with the various eigenmodes of the density 
fluctuations, but not the eigenfrequencies of the system. We have shown in a 
one-dimensional model that the energy- weighted sum rule calculated accord- 
ing to this prescription has exactly the same value as for normal, uncorrelated 
systems, thus we conclude that our prescription eliminates all the spurious 
strength introduced by the constant-A approximation. 

In a simplified model of nuclei, the effects of pairing correlations on the 
isoscalar strength functions has been studied in detail for the quadrupole and 
octupole channels, in the region of giant resonances. In both cases the effects 
of pairing are rather small. More sizable effects are found at lower excitation 
energy, in the region of surface modes which have not been included in the 
present model, but will certainly be more affected by pairing correlations. 



Appendix 

In this Appendix we show that the correlated zero-order response function 
given by the modified density fluctuation (99) satisfies the same energy-weighted 
sum rule (EWSR) as the uncorrelated response function. We assume that par- 
ticles move in one- dimensional square-well potential, so that formulae become 
simpler because the particle velocity does not depend on position: v{e,x) ~ 




Uncorrelated sum rule 



The uncorrelated propagator [7] 




(135) 



2 / deF'ie) J] 



2nujo cos[na;or(a;)] 



1 cos[nujQT{x')] 



n 



T v{e,x) 



cu — nu!o{e) + ie v{e,x') 
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gives the uncorrelated strength function 5''^ (a;) = — -Im / dxdx'Q{x)D'^{x,x\uj)Q{x') 
and the first moment 



oo 





oo 



2 / E (136) 



n>0 



Correlated sum rule 



The modified density fiuctuation (99) allows us to evaluate the correlated 
propagator through the relation 

5q{x, u;)^(3 J dx'D°(x, x', u;)Q(x') , (137) 



givmg 



D^{x,x',uj)= (138) 
-2nu!o /nuJo\cosnu!oT{x) 1 cosnu!oT{x') 



v{e, x) uj — ujn + is v{€, x') 
and the correlated first moment 



n>0 



The only difference between this expression and Eq. (136) is in the form of 
F'(e), which is proportional to a 5-function in (136), while it is smoother in 
the correlated case, however, if the parameter /i is determined by the one- 
dimensional version of Eq. (27), then it can be easly found that, for a square- 
well mean field. 

Ml = Ml . (140) 



The detailed argument goes as follows: both for correlated and uncorrelated 
systems, the number of particles is given by 

oo 

A = J dxdpF{e) = J deT{e)F{e) , (141) 
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with F(e) = 2^^(eF — e) for uncorrclated fermions and F{e) — ^^Po{() in the 
correlated case, while the moments (136, 139) are given by 



j deF'{e)G{e) , 



(142) 



with 



G(e) 



T{e) 



n>0 



(143) 



in both cases, F'{e) obviously differs in the two cases. 
Integrating by parts the last expression in (141), gives 



A^J{e)F{e)\^-jF\e)J{e), 



(144) 



with 



J(e) = j deT{e) 



(145) 



For a square-well potential of size L: 



T(e) = . 
J(e)=- 



/2mL 



Since 



lim J{e)F{e) = hm J{e)F{e) = , 



(146) 
(147) 

(148) 



both for the correlated and uncorrelated distributions, we have 

oo 

A^- J deF\e)J{e). 



(149) 



The explicit expressions of J(e) and T(e) give 
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oo 




for both distributions. From these relations follows that, for a square- well mean 
field, the relation (140) is exact. If we had used the density fluctuation (95), 
instead of (99), to evaluate the correlated propagator, we would have obtained 
a different value of the first moment because the additional term (106) in 
the phase-space density gives an extra contribution to the density response 
function and hence to the EWSR. Because of the fundamental character of 
the EWSR, as well as of the continuity equation, we conclude that this term 
is a spurious contribution generated by the constant-A approximation. 
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